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Abstract 

The longitudinal profile of the most energetic cosmic-ray air shower measured so far, 
the event recorded by the Fly's Eye detector with a reconstructed primary energy 
of about 3.2 • 10 20 eV, is compared to simulated shower profiles. The calculations 
are performed with the CORSIKA code and include primary photons and different 
hadron primaries. For primary photons, preshower formation in the geomagnetic 
field is additionally treated in detail. For primary hadrons, the hadronic interaction 
models QGSJET 01 and SIBYLL 2.1 have been employed. The predicted longitu- 
dinal profiles are compared to the observation. A method for testing the hypothesis 
of a specific primary particle type against the measured profile is described which 
naturally takes shower fluctuations into account. The Fly's Eye event is compatible 
with any assumption of a hadron primary between proton and iron nuclei in both 
interaction models, although differences between QGSJET 01 and SIBYLL 2.1 in 
the predicted profiles of lighter nuclei exist. The primary photon profiles differ from 
the data on a level of ~ 1.5a. Although not favoured by the observation, the primary 
photon hypothesis can not be rejected for this particular event. 



1 Introduction 



Identifying the primary particle type might provide a key to understanding the 
origin of the extreme high-energy cosmic rays (EHECR) with energies around 
and exceeding 10 20 eV. In "bottom-up" acceleration scenarios, hadrons are the 
favoured particles at production site with different predictions ranging from 
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proton to iron-dominated fluxes. On the contrary, in "top-down" scenarios 
generally a large fraction of the observable events are predicted to be photon- 
initiated. In these scenarios, a decay of supermassive "X-particles" is assumed 
which supposedly were created by topological defects such as cosmic strings 
and magnetic monopoles or produced as metastable particles in the early 
Universe. In the decay chain of the "X-particles", a large fraction of EHE 
photons is predicted. For a review, see e.g. [1] and references given therein. 
It should be noted that also in "bottom-up" scenarios a small fraction of 
the particles entering the atmosphere might be photons due to interaction 
processes during the hadron propagation. Also "top-down" scenarios with a 
reduced fraction of predicted EHE photons exist. However, the determination 
of the nature of the EHECR and especially conclusions on the fraction of 
EHE photons will provide strong constraints on models of cosmic-ray origin 
and help to decide between "top-down" and "bottom-up" scenarios. 

Based on an analysis of muons in air showers observed by the Akeno Giant Air 
Shower Array (AGASA), upper limits on the photon flux were estimated to be 
28% above 10 19 eV and 67% above 10 19 - 5 eV at a confidence level of 95% [2]. 
By comparing the rates of near-vertical showers to inclined ones recorded by 
the Haverah Park shower detector, upper limits of 48% above 10 19 eV and 
50% above 4 • 10 19 eV (95% C.L.) were deduced [3]. 

The most energetic cosmic ray reported so far was measured by the Fly's Eye 
detector [4,5] located at Dugway, Utah (40° N, 113° W). The observation was 
performed with a compound eye of 880 photomultiplier tubes monitoring the 
sky. The tubes were arranged at the focal planes of 67 mirrors, each mirror 
having 1.5 m diameter. An air shower is observed by the Fly's Eye as fluo- 
rescent light source where the local light intensity is closely connected to the 
charged particle content of the shower, thus allowing to derive a longitudinal 
shower profile. 

The record event was detected on October 15, 1991, at 7:34:16 UT with a re- 
constructed energy of about (3.2±0.9) • 10 20 eV [5]. The main characteristics of 
this event are summarized in Table 1. Comparing the reconstructed depth of 
shower maximum with expectations from air shower simulations, it has been 
concluded that the event agrees well with hadron-initiated showers. The best 
fits were reached for mid-size nuclei, but due to the experimental uncertain- 
ties and due to shower fluctuations, even nucleons or heavy nuclei were not 
excluded [5]. Similar conclusions were obtained in [6]. For the hypothesis of a 
primary photon, it was pointed out that the event direction is nearly perpen- 
dicular to the local geomagnetic field [5,6,7]. Thus, the possibility of creating 
an electromagnetic cascade ( "preshower" ) in the magnetosphere before enter- 
ing the atmosphere has to be taken into account, with considerable effects on 
the longitudinal shower profile. In a more quantitative estimate, the predicted 
cascade curves for primary photons were claimed to be in disagreement with 
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Table 1 

The main shower parameters of the most energetic Fly's Eye event [5]. The com- 
bined uncertainty is obtained adding the statistical and systematic uncertainties in 
quadrature. The azimuth angle is given counterclockwise from east for the incoming 
direction. 





Best-Fit 


Statistical 


Systematic 


Combined 


Parameter 


Value 


Uncertainty 


Uncertainty 


Uncertainty 


Energy [10 18 eV] 


320 


+35 
-40 


±85 


+92 
-94 


X m ax [g/cm 2 ] 


815 


+45 
-35 


±40 


+60 

-53 


zenith angle [deg] 


43.9 


+1.4 
-0.6 


±1.2 


+1.8 
-1.3 


azimuth angle [deg] 


31.7 


+4.0 
-6.0 


±1.2 


+4.2 
-6.1 



the measurements [6]. 

The preshower calculation in [6] was performed based on [8]. Recently, a nu- 
merical error in [8] has been noted that affects the evaluation of photon con- 
version probability [9]. An independent cross-check of the conclusion of [6], 
that are widely referred to in the literature (see, e.g., [1,10]) seems desirable. 

In this work, the reconstructed longitudinal profile of the Fly's Eye event 
is analyzed employing the CORSIKA shower simulation code [11], which is 
linked to the PRESHOWER program [9] in case of primary photons. The fact 
already expressed in [5] that the shower maximum depth does not uniquely 
identify the type of the primary particle, can not be overcome. However, an 
approach is described here to give an approximate probability of the measured 
event to be consistent with the expectation of a specific primary particle type. 
This approach will naturally take air shower fluctuations into account, which 
in the hadronic case are known to be larger for the lighter nuclei. Two different 
hadronic interaction models, QGSJET 01 [12] and SIBYLL 2.1 [13], are used 
in the comparison. Special attention is paid to the investigation of the primary 
photon hypothesis. 

The plan of the paper is as follows. The calculation of the longitudinal shower 
profiles is described in Section 2. This includes a brief description of the 
preshower physics and of the CORSIKA code. The results are given in Sec- 
tion 3, both for hadrons using different primary masses and interaction models 
and for primary photons. The depth of shower maximum is investigated in Sec- 
tion 3.1. Based on the complete longitudinal profile, the level of (dis-) agree- 
ment between simulations and observation is tentatively quantified in Sec- 
tion 3.2. Conclusions are discussed in Section 4. 
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It should be pointed out that the current paper investigates the possibility of 
assigning a particle type to the event profile published in Ref. [5], taking the 
reconstruction results at face value. No attempt is made to analyze the recon- 
struction method itself; for instance, no conclusion about the reconstructed 
primary energy is made. Also we do not comment on potential astrophysi- 
cal sources of this event as well as particle dependent propagation effects on 
the way from the source to the observer; this has been discussed in detail 
elsewhere [14]. 



2 Calculation of longitudinal profiles 

2. 1 Preshower formation 

At energies above 10 19 eV in the presence of the geomagnetic field, a photon 
can convert into an electron-positron pair before entering the atmosphere. For 
a small path length dr, the conversion probability can be determined as [8] 

Pconv(r) ~a(x(r))dr (1) 
where a(x) = 0.5(a em m e c/h)(B ± /B cr )T(x), X = 0.5(/w/ /m e c 2 )(B L / B 

cr ) j ^errt 

is the fine structure constant, B± is the magnetic field component transverse to 
the direction of photon motion, B cr = m 2 e (? /eh = 4.414- 10 13 G, and T(x) is the 
magnetic pair production function which is negligible if x <S 1, has a maximum 
around x — 5 and then decreases slowly to zero [8,9]. The resultant electrons 
will subsequently lose their energy by magnetic bremsstrahlung (synchrotron 
radiation). The probability of emitting a bremsstrahlung photon by a single 
electron over a small distance dr is given by [15]: 




o 



where I(B±,E,hu) is the spectral distribution of radiated energy, E is the 
electron energy and hv is the energy of the emitted bremsstrahlung photon 
[15,9]. If the energy of the emitted photon is high enough, it can create another 
electron-positron pair. In this way, instead of the primary high-energy pho- 
ton, a number of less energetic particles, mainly photons and a few electrons, 
will enter the atmosphere. We call this cascade a "preshower" since it origi- 
nates and develops above the atmosphere, i.e. before the "ordinary" shower 
development in air. 
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Fig. 1. Energy spectrum of the preshower particles (dotted line) and spectrum 
weighted by energy (solid line) for the conditions of the Fly's Eye event (average of 
1000 simulation runs). 



Preshower features have been investigated by various authors, see for instance 
[6,7,9,16,17,18,19,20,21,22,23]. A detailed description of the code applied in 
the present analysis is given in [9] . In brief, the geomagnetic field components 
are calculated according to the International Geomagnetic Reference Field 
(IGRF) model [24]. The primary photon propagation is started at an altitude 
corresponding to a preshower path length of five Earth's radii. The integrated 
conversion probability at larger distances is sufficiently small. The photon con- 
version probability is calculated according to Eq. (1). For the Fly's Eye event 
conditions, all simulated primary photons converted. After conversion, the re- 
sultant electrons are checked for bremsstrahlung emission with a probability 
distribution of the emitted photon energy following Eq. (2). A cutoff of 10 12 eV 
is applied for the bremsstrahlung photons, as the influence of photons at lower 
energies is negligible for the air shower evolution. The preshower simulation is 
finished when the top of atmosphere is reached. Then, all preshower particles 
are passed to CORSIKA. The resultant energy spectrum of the preshower par- 
ticles is shown in Figure 1. The photons and electrons reach the atmosphere 
with energies below 10 20 eV; most of the initial energy is stored in particles of 
~ 10 19 eV. 
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2.2 Event generation with CORSIKA 



CORSIKA, a standard tool for Monte Carlo shower simulation, has been used 
for calculating the atmospheric shower profiles (version 6.16). In case of pri- 
mary photons, the resultant atmospheric cascade is simulated as a superposi- 
tion of subshowers initiated by the preshower particles. Electromagnetic inter- 
actions are treated by the EGS4 code [25], which has been upgraded to take 
photonuclear reactions as well as the Landau-Pomeranchuk-Migdal (LPM) ef- 
fect [26] into account [11,27]. The LPM effect leads to an increase of the mean 
free path of electromagnetic particles and enhances an asymmetric energy dis- 
tribution of the secondary particles. In particular for electromagnetic particles 
at energies exceeding 10 19 eV, the shower development can be considerably re- 
tarded, and shower-to-shower fluctuations are large. For the primary energy 
and direction of the Fly's Eye event, however, the preshower formation in case 
of primary photons transfers the energy to particles which are less (or not at 
all) affected by the LPM effect. 

To illustrate the importance of accounting for the preshower formation, longi- 
tudinal profiles of primary photons without preshower simulation are displayed 
in Figure 2 together with results of the complete calculation. It is interesting 
to note that the depth of shower maximum in the full simulation is signif- 
icantly smaller also compared to the case of neglecting both preshower and 
LPM effect: Most shower energy is stored in particle energies 1.5 orders of mag- 
nitude below the initial one (see Figure 1). Thus, with an elongation rate of 
~85 g/cm 2 per primary energy decade for unaffected photon showers, a rough 
estimate yields differences of ~130 g/cm 2 in depth of shower maximum, in 
good agreement with the simulation. 

Hadron-initiated showers were simulated with the hadronic interaction mod- 
els QGSJET 01 [12] and SIBYLL 2.1 [13] to study the difference in primary 
type assignment due to different modeling of high-energy hadronic interac- 
tions. Simulation thresholds for the kinetic particle energies of 250 keV for 
electromagnetic particles and 100 MeV for hadrons and muons were adopted. 
The thin sampling option was used which allows Monte Carlo profile calcula- 
tions at highest energies within acceptable computing times. In this option, 
the bulk of secondary particles produced in an interaction is discarded, and 
only "representative" particles are sampled. Weight factors are assigned to the 
selected particles to keep energy conservation [29]. Thin sampling is activated, 
however, only at particle energies much smaller than the primary energy. This 
ensures an undisturbed modeling of the shower profile fluctuations, in par- 
ticular those at the initial states of the shower cascade. Additionally, upper 
limits for the particle weights have been chosen in CORSIKA to keep the con- 
tribution to the profile of an individual particle sufficiently small [30,31,32]. 
More specifically, thinning is activated for particles with energies < 3- lO -6 ^, 
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Fig. 2. Examples of longitudinal shower profiles of primary photons for the Fly's 
Eye event conditions with and without simulating preshower formation in the geo- 
magnetic field. 

with Eq being the primary energy, and particle weights are limited to < 10 6 , 
which is negligible compared to the total number of particles in the shower 
development range relevant for the current analysis. 

Based on the reconstructed values of primary energy and direction quoted 
in Table 1, simulations were performed for primary protons, carbon and iron 
nuclei, both for QGSJET 01 and SIBYLL 2.1, and for primary photons (using 
QGSJET 01 for the small amount of hadron-hadron interactions). For each 
primary particle setting, 1000 events were generated. Additionally, primary 
energy and direction were varied within the given reconstruction errors. For 
all primary particle types it turned out that the calculated profiles showed no 
significant changes, apart from the normalization in case of primary energy 
variations, within the parameter range allowed by the reconstruction uncer- 
tainties. This is expected for hadron-initiated events. Given an elongation rate 
of about 50—60 g/cm 2 per primary energy decade, a 10% change in energy 
corresponds to a shift of shower maximum of ~2 g/cm 2 , which is well within 
the measurement uncertainties and intrinsic shower fluctuations. Primary pho- 
ton profiles are very sensitive to variations of the primary energy or arrival 
direction close to the threshold for preshower formation. The Fly's Eye event 
conditions, however, are well above the preshower onset also when varying the 
primary energy and direction within the reconstruction uncertainties. 
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Fig. 3. Shower maximum distributions of primary hadrons (iron — dashed, carbon — 
dotted, proton — solid line), both for QGSJET 01 and SIBYLL 2.1. The measured 
value is indicated with the combined uncertainty, cf. Table 1. 

3 Comparison of measurement and simulation 



3. 1 Depth of shower maximum X„ 



The distributions of depth of shower maximum for various nuclear primaries 
and for the hadronic interaction models QGSJET 01 and SIBYLL 2.1 are 
shown in Figure 3. The reconstructed depth of the Fly's Eye event of 815 g/ cm 2 
corresponds to a predicted average value of a mid-size nucleus. But due to the 
measurement uncertainties and due to the intrinsic shower fluctuations, any 
of the hadron primaries considered in both models could have produced the 
observed shower depth. 

The average depth of shower maximum and its fluctuations are listed in 
Table 2. While the predictions for heavy nuclei agree for the two models, 
the difference for primary proton amounts to about 35 g/cm 2 . Compared 
to SIBYLL 2.1, composition analyses based on QGSJET 01 would deduce 
somewhat lighter nuclei as primary particles using average shower maximum 
depths. However, for the current investigation this is of minor importance. 
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Table 2 

Average depth and RMS of shower maximum for the given primary particle as- 
sumptions and discrepancy A to the observed shower maximum in units of standard 
deviations. 







QGSJET 01 


SIBYLL 2.1 




photon 


p C Fe 


p C Fe 


X m ax [g/cm 2 ] 


937 


848 808 783 


882 824 785 


RMS{X max ) [g/cm 2 ] 


26 


54 30 22 


47 27 19 
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Fig. 4. Shower maximum distribution of primary photons compared to the recon- 
structed value of the Fly's Eye event. The measured depth is shown with the com- 
bined uncertainty of la c and 2a c , cf. Table 1. 

In Figure 4, the X max distribution of primary photons is displayed. The Fly's 
Eye event depth of maximum of 815 g/cm 2 is smaller than the average X max of 
primary photons of 937 g/cm 2 . Given a combined uncertainty of 60 g/cm 2 for 
the measured X max , the difference between these values is about two standard 
deviations. 



To account also for shower fluctuations, the following method is applied to 
quantify the level of agreement between a primary particle hypothesis and 
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the observation. For each profile % with its individual shower maximum, the 
corresponding probability Pj of consistency with data is determined. The total 
probability P upon which the decision is based whether to reject a specific 
primary particle hypothesis is then calculated from all n s = 1000 profiles 
simulated per primary type as 

P=-£p. (3) 

In this approach, also possible non-Gaussian shower fluctuations predicted 
by the simulations are naturally taken into account. An incorrect considera- 
tion of shower fluctuations might introduce a significant bias in the primary 
type assignment, since the fluctuations are different for the various primary 
particles. 

For primary photons, adopting the combined uncertainty of the data given in 
Table 1, a value of Px max — 6% is obtained. This corresponds to a discrepancy 
of A ~ 1.9cr, where A is defined by P(x 2 > A 2 ) = 6% from the x 2 distribution 
for one degree of freedom. The discrepancy between data and calculation is 
listed in Table 2 for the different primary particle assumptions. 

It might be instructive to illustrate the importance of reducing the measure- 
ment uncertainty. Assuming a resolution in X max of < 30 g/cm 2 as envisaged 
for the detectors of the Pierre Auger Observatory [33,34], the discrepancy of 
primary photons to the reconstructed value would amount to > 3.3<r instead 
of 1.9(7. 

3.2 Complete longitudinal profile 

For the various hadron primaries, differences in the average number of particles 
(shower size) at maximum are known to be on the level of only a few percent 
and thus much smaller than the uncertainty in primary energy reconstruction 
of the Fly's Eye event. In case of unconverted photons, the number of particles 
at maximum could be considerably reduced, which is accompanied by a change 
in the shape of the profile, as can be seen in Figure 2. For the specific conditions 
of the Fly's Eye event, however, all simulated photons converted, and the 
differences in average shower size at maximum compared to primary hadrons 
are small (<8%). Also the profile shapes of the different primaries are in 
reasonable agreement to the data, as shown in Figure 5. Iron-initiated showers 
show slightly narrower profiles than, e.g., proton-induced ones. The difference 
is small, however, compared to the measurement uncertainties. 

It seems worthwhile to mention that a quantity such as "number of particles" 
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Fig. 5. Profile shape of simulated showers initiated by photon, proton, and iron 
primaries (calculated with QGSJET 01) in comparison to the data. Each profile has 
been divided by the individual maximum number of particles and has been 

shifted by the respective depth of shower maximum X % max . 

would require a precise definition if conclusions were attempted based on a di- 
rect comparison of shower sizes provided by different approaches. For instance, 
the number of tracked particles in Monte Carlo calculations depends on the 
simulation energy threshold [35]. Also, shower particles in three-dimensional 
calculations show an angular spread around the shower axis [36,37], which 
should be kept in mind when comparing to an "effective" particle number in 
a one-dimensional calculation. However, due to the relatively large primary 
energy uncertainty of Fly's Eye event, these items are of minor relevance for 
the present analysis. 

The complete longitudinal event profile is displayed together with simulated 
hadron-induced showers in Figure 6 and with photon-initiated events in Fig- 
ure 7. Only a subset of the simulated events, picked at random, is plotted. The 
normalization of the calculated profiles has been chosen such that the number 
of particles at shower maximum roughly coincides with the reconstructed one 
of the Fly's Eye event. More specifically, all CORSIKA shower size profiles 
are scaled by a factor 1.06 in case of primary photons and by 1.10 for primary 
hadrons, which is well within the reconstruction uncertainty of the primary 
energy. 
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Fig. 6. Typical longitudinal profiles of hadron- initiated showers compared to the 
data for different combinations of primary particle and hadronic interaction model 
as assigned. 



The profiles of the different hadron primaries show reasonable agreement with 
the data for both interaction models. In case of primary photons (Figure 7), 
differences between data and simulation are visible that can be attributed to 
a shift in atmospheric depth of the whole profile. Since the measured data 
points are constrained by a common geometry fit, it is important to note that 
their reconstructed atmospheric depth values are correlated [5] . Qualitatively, 
a shift in depth of each point by about 1.5 units of its standard deviation 
would results in a reasonable agreement of data and the photon profiles. In the 
following, a quantitative treatment for testing a primary particle assumption 
against the observed profile is described. 
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Fig. 7. Typical longitudinal profiles of photon-initiated showers compared to the 
data. 

Shower fluctuations are taken into account as discussed in Section 3.1 by av- 
eraging the probabilities Pj for the individual profiles to obtain an overall 
probability P for each considered primary particle type (Eq. 3). The calcula- 
tion of Pi for each individual simulated profile is based on the reconstructed 
profile data as shown in Figure 7 assuming Gaussian probability density func- 
tions for the quoted uncertainties dNj and dXj of the data points. For the 
uncertainties dXj, we take the point-to-point correlation into account and al- 
low all data points to be shifted in the same direction by the same fraction 
/ • dXj when varying the reconstructed geometry. Given a simulated profile i, 
for a certain correlated shift of the rid = 12 data points in horizontal direction 
by fk ■ dXj, a value xfk * s then calculated as 



with ANjk = Nj — N S i m (Xj + fk ■ dXj) denoting the difference between 
the reconstructed event value Nj and the simulated value N sim at depth 
Xj + fk ■ dXj. Factors fk = — 3.0, +3.0 are adopted with a stepsize of 
Af = fk+i — fk = 0.05. From xik-> ^ ne probability P ik that the i-th profile 
fits the shifted data is determined as 

Pik = P(xl, ge , > Xl) (5) 
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Fig. 8. Distribution of probabilities Pi for photon primaries, calculated as described 
in the text. 

by randomly generating a series of artificial data sets according to the uncer- 
tainties dNj around N sim (Xj) and shifting each data set in the same way as 
the real data to obtain the distribution of Xik, g en values [38] . The probability 
Pi is then given by the sum 



Pi=J2 w k- P 



ik 



(6) 



where the individual P^ are weighted for step k according to a Gaussian 
probability density function, 



Wk 



2tt 



J exp(-z 2 /2) dz 



(7) 



fk-h^f 



The resulting distribution of Pj for the photon profiles is plotted in Figure 8. 
Most profiles reach values well above 5%. Following Eq. (3), an average prob- 
ability P ~ 13% is obtained. This corresponds to a discrepancy between pho- 
tons and data of about 1.5a. The results obtained for the primary photon and 
hadron hypotheses are summarized in Table 3. 
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Table 3 

Probability P of a given primary particle hypothesis to be consistent with the ob- 
served Fly's Eye event profile and corresponding discrepancy A in units of standard 
deviations. 



QGSJET 01 SIBYLL 2.1 





photon 


P 


C 


Fe 


P 


c 


Fe 


P[%] 


13 


43 


54 


53 


31 


52 


54 


A [a] 


1.5 


0.8 


0.6 


0.6 


1.0 


0.6 


0.6 



4 Discussion 



For any considered primary hadron assumption and for both hadronic inter- 
action models, the values given in Table 3 confirm that no hypothesis can 
be rejected, as expected from the X max investigation and in agreement with 
previous considerations in the literature [5,6]. Also in the framework of the 
SIBYLL 2.1 model, primary protons, which on average reach the maximum 
~70 g/cm 2 deeper in the atmosphere than the observed profile, are not ex- 
cluded due to the measurement uncertainties and shower fluctuations. They 
are only slightly less favoured (discrepancy of 1.0a to data) with respect to 
heavier nuclei (0.6er) or primary protons in the QGSJET 01 model (0.8a). 

Compared to primary hadrons, photons are less favoured. With a discrepancy 
to the Fly's Eye event profile on the level of 1.5<r, however, also the primary 
photons are not excluded by the data. 

This result does not confirm a previous analysis that claimed the primary 
photon hypothesis to be inconsistent with the experimental observations [6]. 
Compared to our calculations, the shower maximum of primary photons had 
been obtained at depths larger by ~90 g/cm 2 . Part of this difference might 
arise from the preshower calculation, as might be expected due to the numer- 
ical error in Ref. [8] pointed out in Ref. [9]. Our analysis yields an average 
energy fraction carried by particles with energies above 10 19 eV of ~ 40% 
(compare Figure 1). As this is smaller than the value > 60% given in [6], a 
"faster" shower development is expected, also due to the LPM effect being 
less important for smaller particle energies. 

The discrepancy between photons and data might be diminished below 1.5a 
by additional uncertainties not yet considered in the quantitative treatment. 
For the simulated profiles, no "theoretical" uncertainty has been assumed 
when testing the primary photon hypothesis. Though much less affected by 
high-energy extrapolations of hadronic interactions than hadron primaries, un- 



15 



certainties are associated, for instance, with the photonuclear cross-section at 
the highest energies. Given the specific conditions of the Fly's Eye event, pho- 
tonuclear reactions seem of minor importance. For photons around 10 19 eV, 
the expected photonuclear cross-section [39] is more than two orders of mag- 
nitude below the pair production cross-section. As preshowering additionally 
distributes the initial primary energy among many particles entering the at- 
mosphere, photonuclear interactions of a small subset of these particles (of 
order 1 out of 300 particles) impose only a marginal effect on the overall 
shower profile for the simulated primary photons. Allowing, however, for a 
larger increase of the photonuclear cross-section with energy, photon-initiated 
events would become more hadron-like. The expected photon profiles would 
be shifted towards the observed one which would reduce the discrepancy. 

For the reconstructed data, recent studies of the atmospheric density pro- 
file [40] indicate a significant additional uncertainty in the reconstructed at- 
mospheric depth. The air density profile is required for converting the altitude 
of a shower track observed by the telescope to its derived atmospheric depth, 
as the physics interpretation is mostly based on the latter quantity. Without 
proper correction, variations of the density profile might lead to misinter- 
pretations of the reconstructed atmospheric depth of 10-20 g/cm 2 or more, 
depending on the specific conditions of the observation. Such density profile 
variations have been confirmed in investigations of atmospheric data taken 
at Salt Lake City airport, i.e. not far from the Fly's Eye site, that are avail- 
able starting from 1996 [41]. Moreover, a study of October density profiles 
measured over several years at Salt Lake City airport indicates a systematic 
underestimation of about 10—15 g/cm 2 on average in X max when using the 
US standard atmosphere for geometries similar to the Fly's Eye event [42]. 
Assuming such a shift in atmospheric depth of the reconstructed profile, the 
discrepancy between photons and data would be reduced by ~0.2er. 

Even without assuming a re-analysis of the Fly's Eye event that would push 
the data closer to the photons, it seems fair to conclude that additional uncer- 
tainties exist so that the potential to exclude a specific primary particle type 
hypothesis is reduced. This strengthens the conclusion that photons are not 
ruled out as primary particle for the Fly's Eye event. 

Clearly, a larger statistics of high-energy events is necessary to allow strin- 
gent constraints on the photon hypothesis. Let us illustrate the sensitivity of 
the Pierre Auger Observatory [33,34] for giving an upper limit on the photon 
fraction in case no photon detection can be claimed. Depending on the actual 
high-energy particle flux and including a duty cycle of 10—15%, the fluores- 
cence telescopes of the Auger experiment are expected to record about 30—50 
shower profiles with primary energies exceeding 10 20 eV within a few years of 
data taking. For simplicity, let us consider a constant probability e for each of 
these n observed profiles to be originated by a photon. Then, the probability 
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Table 4 

Simulated average muon number and RMS on observation level for the Fly's Eye 
event conditions. 







QGSJET 01 


SIBYLL 2.1 




photon 


p C Fe 


p C Fe 




3.1 


12.9 14.0 15.3 


9.6 11.0 11.8 


RMS(iV M /10 8 ) 


0.3 


1.5 0.8 0.6 


1.2 0.6 0.5 



that n 7 out of n profiles were photon-initiated is 

%) = 6^(l-eP( n ;) . (8) 

For a choice of n = 40 and e = 5%, with 95% confidence level a photon 
fraction exceeding 10% could be excluded, which would provide an important 
constraint for explaining the EHECR origin. In case of the Fly's Eye event, 
a value of e = 5% corresponds to reducing e.g. the uncertainties dXj by a 
factor 1.5, which seems well in reach for the Auger experiment. The sensi- 
tivity level from the profile measurements alone might even be increased, as 
the most energetic events will be observed from different telescope positions. 
Further improvement might be gained by including data from the Auger sur- 
face detectors. Apart from the virtually 100% duty cycle of this surface array, 
information on the muon content of the shower would additionally strongly 
constrain potential primary photons. The average muon number at observa- 
tion level for the events simulated in the current analysis is given in Table 4. 
In case of photon primaries, the muon number is reduced by more than a fac- 
tor 3 relative to hadron primaries. It is also interesting to note that the muon 
numbers predicted by QGSJET 01 and SIBYLL 2.1 differ quite significantly 
from each other. Compared to the overall profile, the muon number is much 
more sensitive to different extrapolations of hadronic interaction features to 
the highest energies. Combining the various observables will help to improve 
the modeling of hadron-initiated showers. 
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